Nanoscale Damage During Fracture in Silica Glass 
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Abstract. We report here atomic force microscopy experiments designed to uncover the nature of 
failure mechanisms occuring within the process zone at the tip of a crack propagating into a silica 
glass specimen under stress corrosion. The crack propagates through the growth and coalescence of 
nanoscale damage spots. This cavitation process is shown to be the key mechanism responsible for 
damage spreading within the process zone. The possible origin of the nucleation of cavities, as well as 
the implications on the selection of both the cavity size at coalescence and the process zone extension 
are finally discussed. 
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1. Introduction 

Slicate glasses are often considered as the archetype of brittle materials. It is commonly 
thought (Kelly et al, 1967; Lawn et al, 1980; Guin and Wiederhorn, 2005) that its 
fracture is similar to cleavage - with a crack progressing through sequential bond 
ruptures without involving any damage ahead of the crack tip. However, some recent 
observations calls into questions this scenario: 

— The morphology of fracture surfaces in glasses exhibits scaling features similar to 
the ones observed in a wide range of quasi-brittle and ductile materials, e.g. oxide 
glass, polymers, metallic alloys, wood, rocks... (see Bouchaud (1997) and refer- 
ences therein). This strongly suggests the existence of some underlying generic 
mechanisms within the process zone common to all these materials. 

— The deformation field was shown not to fit with the linear elastic predictions over 
a fairly large region (of the order of a hundred of nanometers) in the vicinity of 
the crack tip (Guilloteau et al., 1996; Henaux and Creuzet, 2000). 
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— Molecular Dynamic (MD) simulations indicate that fracture in silica glass pro- 
ceeds through the growth and coalescence of nanoscale cavities (Van Brutzel et 
al, 2002; Rountree et al, 2002; Kalia et al, 2003). 

These various results have led us to investigate experimentally the failure mecha- 
nisms occuring in glass at its microstructure scale, the nanoscale. The first series of 
experiments, reported in Celarie et al. (2003a, 2003b) and Marliere et al. (2003), were 
carried out on aluminosilicate vitroceramics. They clearly reveal that crack progresses 
through the growth and coalescence of nanoscale damage cavities. We describe here 
our recent studies performed on a minimal elastic vitreous medium, pure silica. The 
experimental setup is described in section 2. In section 3, we present the experimental 
results. As in the aluminosilicate vitroceramics, crack propagation proceeds through 
the nucleation, growth and coalescence of damage cavities (section 3.1). The rate of 
cavity nucleation and the size of cavities at coalescence are shown to set the mean crack 
growth velocity as measured at the continuous scale (section 3.2). This nanocavitation 
process is shown to set the process zone size (section 3.3). Finally, section 4 is devoted 
to a discussion on a possible scenario explaining the existence of damage cavities ahead 
of the crack tip, and its implications as regard to the cavity size at coalescence (section 
4.1) and the process zone size(section 4.2). 



2. Experimental Setup 

The experimental set-up has been described in detail in Prades et al. (2005) and is 
briefly recalled below. Fracture was performed on double cleavage drilled compression 
parallelepipedic samples (size 5 x 5 x 25 mm 3 ) with a cylindrical hole drilled in the 
center (radius 0.5 mm). A gradually increasing uniaxial compressive load was applied 
to the sample (Figure 1). Once the two cracks are initiated symmetrically from the 
hole, the load is held constant. In this geometry, the stress intensity factor Kj can be 
related to the crack length using the expression given by He et al. (1995). While the 
cracks length increases, Kj decreases. Under vacuum, the crack would stop when Kj 
gets smaller than the material fracture toughness Kj c . But, in our room atmosphere 
(relative humidity 45%, temperature 26°C), the corrosive action of water on glass 
allowed slow, sub-critical, crack propagation. The crack tip advance was then slow 
enough to be probed by AFM. 

The crack growth velocities v considered in this study range between 10~ 12 m/s 
and 10~ 6 m/s. In this velocity range, v increases exponentially with Kj (Figure 2) in 
agreement with stress enhanced activated process models as proposed by Wiederhorn 
(1967) and Wiederhorn and Boltz (1970). The velocity could then be tuned by ad- 
justing the external applied load for a measured crack length. It was varied over three 
decades, from 10~ 9 and 10 -12 m/s (the maximum reachable velocity being set by the 
recording time of an AFM frame, around 3 minutes). 

In all the following, the reference frame (e x ,e y ,e z ) is chosen so that e x , e y and 
e z are parallel to the crack propagation, tension loading at the crack tip and sample 
thickness directions respectively. 
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Figure 1. Experimental setup showing the loading DCDC configuration used to fracture the glass 
specimen under stress corrosion. 
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Figure 2. Variation of the mean crack growth velocity v (in m.s -1 ) as a function of the stress intensity 
factor Ki (in MPa.m 1/2 ). The axes are semi-logarithmic. Diamonds and dots correspond to optical 
and AFM measurements respectively. The line is a fit using v = vo exp(K] / Kq ) as expected for stress 
enhanced activated process models. 
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3. Experimental results: Damage spreading mechanisms within the 

process zone 

3.1. Evidence of damage cavities ahead of the crack tip 

Our experimental setup has allowed to observe the crack progression within the process 
zone in real time at the nanoscale. Figures 3(a-f) show six successive AFM topographic 
frames in the vicinity of the crack tip in the pure silica glass. For this specific sequence, 
the crack propagates at an average velocity v ~ 4.10~ n m/s. One can clearly see a 
depression ahead of the crack tip. This cavity grows to a typical size of 100 nm 
in length and 20 nm in width, and then merges with the advancing main crack to 
make it cross the whole area of observation. In other words, the crack front does 
not propagate regularly as commonly stated (Kelly et al., 1967; Lawn et al, 1980; 
Guin and Wiederhorn, 2005), but progresses through the growth and coalescence of 
cavities. This scenario is fully consistent with what was observed both experimentally 
in aluminosilicate vitroceramic under stress corrosion by Celarie et al. (2003a, 2003b) 
and Marliere et al. (2003) (Figure 4(bottom)) and numerically in dynamically breaking 
samples of amorphous Silica by Van Brutzel et al. (2002), Rountree et al. (2002) and 
Kalia et al. (2003) (Figure 4(top)). 

Fracture Surface Topography Analysis (Kobayashi and Shockey, 1987; Miyamoto 
et al., 1990) was then performed to ensure that the spots observed ahead the crack tip 
correspond actually to damage cavities. In a ductile scenario, the growth of damage 
cavities is expected to induce irreversible plastic deformations that will leave visible 
prints on the post-mortem fracture surfaces. We have thus determined the fracture lines 
after the crack has crossed the area of interest (white lines in Fig. 3f ) and reconstituted 
virtually the intact material by placing the line on the right - line R - to the left of 
the line on the left - line L (Bonamy et al., 2005). By translating gradually the line R 
to the right, one reproduced qualitatively the chronology of the cavity growth (Fig. 
3a'-f). This provides a rather strong argument to relate these nano-scale spots to 
damage cavities. 

It is worth to note that the reconstructed cavity is found to be thinner than the 
real one. This can be understood as follows: The shape of the real cavity is given not 
only by the irreversible - permanent - part of the opening displacements at the free 
surface, but also by the reversible - elastic deformation -. This latter vanishes once the 
crack has crossed the area of interest and the stresses have relaxed, - and therefore 
can not be taken into account in the reconstructed frames (Bonamy et al, 2005). In 
others words the remnants of these nanospots growing ahead of the crack tip on the 
post-mortem fracture surfaces are expected to be much smaller - and thus much harder 
to detect - than the spots themselves. This may explain why Guin and Wiederhorn 
(2005) did not succeed to see evidences of the pores' remnants from the analysis of 
the mismatch between the two opposite post-mortem fracture surfaces. 

3.2. Kinematics of damage cavities 

From the sequences represented in Figure l(a-f), the temporal evolutions of the main 
crack tip (CT), the forward front tip (FF) and the backward front tip (BF) of the 
cavity were then determined, and the corresponding velocities extracted (Figure 5). 
These velocities were found to be v CT = 4.10 -12 m/s, v FF = 1.1. 10 -11 m/s and 
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(a') (b') (c') (d') (e') (f) 

Figure 3. (a)-(f) Sequence of topographic AFM frames in the vicinity of the crack tip showing the 
propagation of the crack in stress corrosion regime (v = 4.1CP 11 m.s -1 ). The scan size is 470 x 135 nm 2 
and the height range over 3 nm. The 6 frames were taken during the five hours necessary for the main 
crack to cross the area of interest. (a,b) Appearance of a nanometric damage cavity ahead of the crack 
tip, (b-d) growth of the cavity prior to crack propagation and (f,g) coalescence of the cavity with the 
main crack. The points A, B and C locate the position of the main crack tip (CT), the backward front 
(BF) and the forward front (FF) tips of the cavity respectively, (a')-(f') reconstructed frames using 
the FRASTA method 

v BF = 1.2. 10 -11 m/s for the main crack tip, the forward front tip and the backward 
front tip of the cavity, respectively, i.e. significantly smaller than the mean crack tip 
velocity v = 4.10 -11 m/s as measured optically (Prades et al, 2005). In other words, 
the crack growth velocity as observed at the continuum scale is dominated by the 
accelerating phases corresponding to cavity coalescence with the main crack front. 

3.3. Nanocavitation and process zone 

We then determined the extension of the process zone (Bonamy et al, 2005; Prades 
et al, 2005). Linear Elastic Fracture Mechanics (LEFM) predicts that the stress field 
exhibits a universal square root singularity at the crack tip in an intact linear elastic 
medium. One thus expects a square root singularity in the components of the stress 
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Figure 4- Top: Snapshots of atoms as observed in a MD simulation of dynamic fracture in amorphous 
Silica (from Kalia et al, 2003). Bottom: Sequences of three sucessive AFM snapshots showing the 
vicinity of the crack tip at the surface of an aluminosilicate glass specimen (from Celarie et al., 2003). 
In both cases, the crack progresses through the growth and coalescence of damage cavities 
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Figure 5. (a) Temporal evolution of the main crack front (CF) and the backward front (BF) of 
the cavity and the forward front (FF) of the cavity represented by points A, B and C on the 
frames of Fig. 3(a)-(e) respectively. The velocities of these fronts are determined through linear fits: 



v CF ~ 4.10~ lz m.s -1 , v ar ~ 1.2.10" 11 m.s" 
the velocity at the continuum scale: < v >~ 
at the scale of the process zone. 
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4.10 m.s . (b) Sketch of the crack tip propagation 
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tensors, and in the out-of-plane displacement at the free surface of the specimen. 
The surface topography was thus imaged in the vicinity of the crack tip (Figure 
6a), and compared to the LEFM predictions (Figure 6b). Far from the crack tip, 
one recovers the LEFM predictions. But below a given threshold, the experimental 
curves depart from the predictions. This departure point sets the extent of the process 
zone (Guilloteau et al, 1996; Henaux and Creuzet, 2000; Celarie et al, 2003a). This 
extension was found to be around 100 — 300 nm in length (x direction) and around 
20 nm in width (y direction) in the stress corrosion experiment (y ranging from 8.10 -12 
to 10" 9 m/s) (Bonamy et al, 2005; Prades et al, 2005). It was found to be significantly 
smaller, around 10 — 25 nm in length, and 7 — 11 nm in width for the MD simulation 
of dynamic crack propagation (v ranging from 10 to 400 m/s) (Prades et al, 2005). 
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Figure 6. Measurement of the extension of process zone, (a) 1 x 1 /im AFM topographical frame of 
the vicinity of the crack tip. The crack growth velocity was v ~ 3.10 -11 m/s (b) Variation of the 
out-of-plane displacement u z as a function of the distance r from the crack tip measured along the 
direction of crack propagation (s-axis). The axes are logarithmic. The straight red line corresponds to 
the LEFM predictions. For r smaller than a given value R c ~ 200 nm, the experimental curve departs 
from the LEFM predictions. This value sets the extent of the process zone along the x direction for 
this specific experiment. 



We finally compared the extension of the process zone as defined above to the 
extent of the zone made "porous" in the vicinity of the crack due to the presence 
of cavities. Those were shown to coincide (Prades et al, 2005). This indicates that 
nucleation and growth of cavities provide the dominant mechanism responsible for 
damage spreading within the process zone. 



4. Discussion: On the relevant lengthscales 



Observations of damage cavities in amorphous silica under stress corrosion are con- 
sistent with what was observed experimentally in aluminosilicate vitroceramics under 
stress corrosion by Celarie et al. (2003a, 2003b) and Marliere et al. (2003), and nu- 
merically during dynamic failure of amorphous Silica by Van Brutzel et al. (2002), 
Rountree et al. (2002), and Kalia et al. (2003). This indicates that the origin of such a 
nanoscale damage mode is inherent to the amorphous structure and does not depend 
on the precise glass composition. Origin of the damage cavities should be found in 
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the amorphous structure, which contains inherent density and/or stress fluctuations 
at the nanoscale. 

4.1. Cavity size at coalescence 

In this scenario, the typical size of the damage cavities can be understood as follows. 
Intrinsic amorphous fluctuations induce low toughness regions with a density po of 
weak bonds, i.e. bonds that can be broken ahead of the crack tip. These " nano-cracks" 
behave as stress concentrators and can grow under the stress imposed by the presence 
of the main crack, to give birth to damage cavities. In dynamic fracture, the stress 
within the process zone is sufficient to activate any of these " low toughness sites" . The 
typical distance between two nucleation sites - and consequently the typical size of 
cavities at coalescence - is given by the lengthscale 5q of the fluctuations induced by the 
amorphous structure. Since 5o is of the order of the nanometer, one expects a cavity 
size of the order of one nanometer which is compatible with the MD observations (Van 
Brutzel et ai, 2002; Rountree et ai, 2002; Kalia et al, 2003). 

This scenario can be extended to the stress corrosion regime (Prades et al, 2005). 
In that case, the stress level in the process zone is too low to make low toughness 
sites give birth to damage cavities. However, the corrosive action of water activates 
some of them. The probability p/po - where p refers to the density of activated sites - 
is expected to scale as p/po oc exp(aG//cT) where a G, k and T refer to a typi- 
cal activation area, the mechanical energy within the process zone, the Boltzmann 
constant, and the temperature respectively. Since the crack growth speed v scales 
also as v oc exp(aG/kT) (Wiederhorn, 1967; Wiederhorn and Boltz, 1970; Prades 
et al, 2005), one gets p/po oc v. The typical distance S between two nucleation 
points - and consequently the typical size of the cavity at coalescence - scales as 
5 oc p -1 / 3 . Hence 5/5q oc f -1 / 3 . In glasses under stress corrosion, crack growth stops 
being dominated by the rate of hydrolysis chemical reactions when v is greater than a 
threshold v* classically described as the limit between stress corrosion regions I and II 
in the literature (Wiederhorn, 1967; Wiederhorn and Boltz, 1970; Lawn, 1993). The 
corresponding mechanical energy G(v = v*) is thus expected to be sufficient to activate 
all the low toughness sites: 5(v = v*) ~ do- Finally, one gets 5/5$ ~ (v/v*)^ 1 ^ 3 in the 
ultra-slow stress corrosion regimes. For our experimental conditions (relative humidity 
45%, temperature 26°C), v* ~ 10~ 5 m/s (Wiederhorn, 1967; Wiederhorn and Boltz, 
1970; Lawn, 1993). Since do ~ 1 nm, one gets 5 ~ 100 nm for v ~ 4.10~ n m/s, which 
is in good agreement with observations (Figure 3). 

4.2. Size of the process zone 

Let us now discuss the typical size R c of the process zone. Analysis performed by 
Swiler et al. (1995) suggests the presence of voids as large as 0.45 nm in unstressed 
silica glass. In first approximation, the process zone extent can be assimilated to the 
maximum distance from the crack tip where the stress imposed by the presence of the 
main crack is sufficient to make such a void grow. At this distance, the average stress 
level is expected to be Kj / '\J2itR c (Inglis, 1913). This stress is then concentrated at 
the notches of the voids and becomes a vmds ~ (1 + 2r)Kj / ^2ttR c where r refers to 
the typical aspect ratio of the void. The void will grow when o~ vmds becomes larger 
than the intrinsic strength of the material a*. Hence, R c ~ (l/27r)(l + 2r) 2 (Kj/a*) 2 . 
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The value of intrinsic strength of silica glass measured in vacuum at room tem- 
perature is estimated around a* ~ 10 — 12 GPa (Kurkjian et al., 2003). Voids in 
unstressed silica glass are expected to be around 0.2 — 0.45 nm thick and 0.45 nm 
large (Swiler et al., 1995), which leads to an aspect ratio r ~ 1 — 2. As dynamic 
fracture occurs for Kj = Kj c = 1 MPa.m 1 / 2 , one thus expects a process zone size 
R c ~ 10 — 40 nm, which is in good agreement with the measurements performed on 
the MD simulations (section 3.3). 

The intrinsic strength of Silica glass in humid air was found to be significantly 
smaller, around a* ~ 3 — 4 GPa (Proctor et al, 1967). Since the stress corrosion 
experiments were performed with a stress intensity factor Kj ~ 0.4 — 0.5, one expects 
a process zone size R c ~ 15 — 110 nm. This order of magnitude is compatible with what 
was measured from the AFM topography (Figure 6). One may however wonder how 
water molecules can travel from the crack tip along such a large distance within the 
glass. For water to travel over this distance in typically an hour, it requires a diffusion 
coefficient D ~ 10 -14 cm 2 s _1 (Bonamy et al., 2005) fives orders of magnitude larger 
than what is commonly observed in the absence of stress (D ~ 10~ 19 cm 2 s _1 ). 
However, Tomozawa and Han (Tomozawa et al., 1991) have shown that under stress, 
diffusion is considerably accelerated, and hence it is perfectly plausible that water 
molecules react with Si-0 bonds at a distance ~100 nm from the crack tip within the 
bulk, under our experimental conditions. 

5. Conclusion 

Ultra-slow crack propagation was observed in real time at the nanoscale through AFM 
in amorphous Silica. The main results from our observations are: 

(i) Within the process zone, the crack tip does not propagate regularly, but through 
the growth and coalescence of damage cavities. 

(ii) Velocity of the main crack tip and the cavity tips as measured at the damage zone 
scale are shown to be significantly smaller than the mean crack growth velocity 
as measured at the continuous scale. In other words, crack growth is intermittent 
and dominated by the accelerating phases corresponding to the cavity coalescence 
with the main crack front. 

(iii) The nanocavitation process provides the dominant mechanism responsible for 
damage spreading within the process zone. 

Observations of damage cavities in amorphous Silica under stress corrosion are 
consistent to what was observed experimentally in Aluminosilicate glasses under stress 
corrosion Celarie et al. (2003a, 2003b) and Marliere et al. (2003), and numerically 
during dynamic failure of amorphous Silica (Van Brutzel et al., 2002; Rountree et 
al, 2002; Kalia et al, 2003). This indicates that the existence of this nanoductile 
mode is inherent to the amorphous structure and does not depend on the precise glass 
composition. 

The origin of cavitation should be found in the intrinsic toughness fluctuations 
induced by the amorphous structure of the material. This results in low toughness 
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region sites and/or low density regions that behave as stress concentrators and give 
birth to cavities ahead of the main crack tip. Such a scenario was shown to capture 
quantitatively the cavity size at coalescence and the process zone extension in both 
the dynamic fracture simulations and the stress corrosion experiments. 

This failure mechanism through growth and coalescence of damage cavities is very 
similar to what is observed classically - albeit at other length scales - in a wide range 
of materials, e.g. aluminosilicate vitroceramics (Celarie et al, 2003a, 2003b; Marliere 
et al, 2003), nanophase materials (Rountree et al., 2002; Kalia et al, 2003), PMMA 
(Ravi-Chandar and Yang, 1997) and polymers (Lapique et al, 2002). We argue that 
such mechanism is generic to crack propagation - the main difference resides in the 
typical lengthscales over which cavities are observed. These lengthscales are controled 
by the typical size of the microstructure eventually modified through environmental as- 
sisted activated process like stress corrosion or fatigue. The aspect ratio of the cavities 
- much larger in glasses and quasi-brittle materials than in metallic alloys - reflects 
the ability of the material to deform irreversibly. Such a scenario may explain the 
puzzling similarities observed in the scaling properties exhibited by fractures surfaces 
in a wide range of materials (Bouchaud, 1997). 

Let us finally add that the present study sheds light on the role of the spatial 
fluctuations induced by the material microstructure, but passes over the temporal 
fluctuations in silence. Interaction of a growing crack with the material microstructure 
results in the release of acoustic waves (Ravi-Chandar and Knauss, 1984; Bonamy and 
Ravi-Chandar, 2003) which play a significant role in the energy dissipation properties 
within the process zone. Their understanding represents interesting challenges for 
future investigations. 
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